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Abstract 



The intrinsic pinning in high-temperature superconductor with long range 
antiferromagnetic order of rare-earth ions confined to the isolating planes is 
described. The interaction of antiferromagnetic and superconducting subsys- 
tems in the mixed state may lead to the creation of spin-flop domain along 
the vortex core. It is shown that the behavior of several physical quanti- 
ties such as magnetization, activation energy, current-voltage characteristic 
and flux creep changes when the direction of the external field changes in 
the basal a-b plane. It is also shown that the decay of the trapped flux is 
logarithmic function of time. 



1 INTRODUCTION 



The problem of interplay between long-range magnetism and superconduc- 
tivity was studied thoroughly several years ago in Chevrel phases, rhodium 
borides and rhodium stannides [|l|] . On the contrary, this phenomenon has not 
attracted much attention in the high-temperature superconductors (HTS). 
Both competitive phenomena coexist in classical superconductors because 4f 
electrons of rare-earth (RE) atoms responsible for magnetism and 4d elec- 
trons of molybdenum chalcogenide clusters (rhodium boride or rhodium stan- 
nide clusters) responsible for superconductivity are spatially separated from 
each other. The situation seems to be very similar in layered HTS. Here mag- 
netic order is produced by the regular lattice of RE ions occupying isolating 
layers electrically isolated from the superconducting Cu-0 planes. Therefore 
spin interaction between the local magnetic moments and the conduction 
electrons is to weak to inhibit superconductivity. The typical example of 
that layered system is ErB^CusOy. This compound has tetragonal unit 
cell with small orthorhombic distortion in the a-b plane ||. The Er ions 
form two sublattice antiferromagnetic structure of magnetic moments lying 
parallel and antiparallel to the b direction in the a-b plane ||. Recently 
discovered RE nickel boride-carbides M may serve as an another example of 
layered magnetic superconductors. The structure of these compounds is sim- 
ilar to that of HTS and consists of RE-C layers separated by M2B2 sheets. 
For example in ErNi 2 B 2 C |J the antiferromagnetic structure is associated 
with magnetic moments of Er +3 ions, which order below 6 K in a trans- 
versely polarized planar sinusoidal structure propagating along a or b axis 
with Er moments parallel to the a or b axis respectively. 

In this paper we consider the structure shown on Fig.l that we believe 
simulates real structure of many antiferromagnetic layered superconductors. 
It consists of superconducting layers of thickness d s and magnetic moments 
of RE ions running parallel and antiparallel to the b axis in the isolating 
layers of thickness di, d = d s + di ~ di. 

An antiferromagnet with two sublattices shows different magnetization 
behavior in an external magnetic field applied parallel or perpendicular to 
the easy axis ||. The main difference is that antiferromagnetic configuration 
(AF) is unstable in perpendicular field. When the external field is applied 
parallel to the magnetic moments the AF phase is stable up to the critical 
field Ht above which spin-flop phase develops. This was the basic idea put 
forward to explain strange behavior of DyMoeSs in applied magnetic field J7J. 
Neutron scattering experiments || have shown that this compound develops 
long range AF state in the presence of superconductivity. Unexpectedly 
some ferromagnetic peaks were observed in the mixed state when the applied 
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field exceeded 200 Oe, considerably below Hc2- The idea of induced spin- 
flop transition near the core of the vortex was introduced to explain these 
observations. The ferromagnetic-like order confined along the vortices was 
belived to be responsible for additional peaks above 200 Oe. 

In the present model (see also ||) we come across two following situations: 

1. The external magnetic field is pointing in the b direction and exceeds 
Ht in the vortex core. Then the magnetic subsystem in the core undergoes 
the spin-flop transition. 

2. The external magnetic field is applied along the a direction. Then the 
AF alignment is unstable against spin-flop transformation and the vortex 
state is created in the uniform spin-flop phase. 

Assume the first case. The field intensity in the vortex core in the low 



vortex density regime doubles the value of the external field intensity fLO 
The vortices are created in the uniform antiferromagnetic medium if the 
condition Hc\ < H < \Ht is fulfilled. When the external field reaches \Ht, 
the magnetic field inside the core is Ht whereas outside it decreases from this 
value. In this way the spin-flop domain appears along the vortex. Finally, 
depending on the external field intensity, there can be two types of vortices 
in the specimen. On the other hand, in the case 2, the vortices do not change 
their structure and are similar to that in nonmagnetic HTS. 

The magnetic field aligned parallel to the conducting planes makes the 
vortex lattice to accommodate itself to the layer structure so that the vortex 
cores lie in between the superconducting sheets. A current density j flowing 
along the planes exerts a Lorentz force on the vortices in the c direction. 
Intrinsic pinning barriers are formed on strongly superconducting layers. The 
vortices loose large amount of condensation energy when they move across 
the barriers. Depending on the external field intensity and its direction in 
the a-b plane we deal with vortices of two types. It is interesting therefore to 
see how changes the intrinsic pinning barrier during the creation of magnetic 
domain along the vortex. 



2 LONDON EQUATIONS 

The electrodynamics of layered superconductors is based on the Lawrence- 
Doniach model that in the London approximation gives the following func- 



tional of the free energy of superconducting subsystem [11 




2vr 2 2 



( V $ n + — A p ) 2 + - [1 - cos( X J] d 2 r } (1) 
V?o r j ) 



2 



where Tj = d(A c /A a b),A = (A p , A z ), A p = (A x ,A y ), <p denotes the flux 
quantum, /i the magnetic permeability of the vacuum, A c ,£ c and \ a b,£, a b 
the magnetic field penetration depth and the coherence length parallel and 
perpendicular to the layers respectively. The parameter rj plays the role 
of the effective radius of the vortex core. The term with cosine gives the 
Josephson current due to the gauge invariant phase difference between the 
layers 

X„ = $„+i-$ n + — / A z dz. (2) 

The antiferromagnetic two sublattice subsystem with single ion anisotropy is 
described with the following free energy functional: 

F m = diJ2 f ( JMi„M 2n -KY, M in± - 7S £ (VM„) 2 )/ v+ i(^ 2 r, 

n J { i=l i=l a=x,y,z ) 

(3) 

where M n = Mi„+M 2 „ is the sum of the magnetization vectors of the sublat- 
tices in the n -th insulating layer, M in ± is the component along the anisotropy 
axis of the magnetization sublattice vector in the n-th layer, J denotes the 
exchange constant between two sublattices, K is the single ion anisotropy 
constant, ^7 is the magnetic stiffness length, the factor f n ,n+i(z) = 1 for 
insulating layers and otherwise, and M = |M ln | = |M 2n |. Finally we add 
the magnetic field energy to obtain the free energy of the entire system 

F = F s + F m + ^-J H 2 d 2 rdz. (4) 

The coupling between the magnetic and superconducting subsystem is sup- 
posed to be electromagnetic. This means that both order parameters are 
coupled through the vector potential A 

B = rot A = yU H + M, (5) 

|^| 2 wj, = |rotH|, (6) 

where B ( in MKSA units) is the vector of a magnetic flux density ( magnetic 
induction) and H is the vector of a thermodynamic magnetic field intensity. 
The analysis of a stability of the functional (3) against the spin-flop tran- 
sition gives the thermodynamical critical field of the upper limit of the AF 
configuration equal to 



H T = M q JK(J — K). 
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The London equations resulting from (4) are following: 



B x + A 2 .— rot^B - M) - A 2 o6 — rot,(B - M) = <p 6(y)6(z), 

B y + A 2 6 — rot^B - M) - A c 2 — rot,(B - M) = 0, (7) 

B z + A 2 6 £rot,(B - M) - A 2 b |-rot x (B - M) = 0. 

This equations should be supplemented by the appropriate set of differential 
equations describing the spatial distribution of magnetization. Simple con- 
jecture, as shall be seen later, can make the calculations less complex and 
at the same time does not oversimplify the problem. The London model 
of continuous superconductors may be used at length scales larger than the 
coherence length, i.e. the core dimension. The structure of a vortex lying in 
the a-b plane in a layered superconductor with Josephson coupling between 
adjacent layers resembles the Abrikosov's one except that the order param- 
eter does not vanish anywhere [[EJ. Instead there exists a region, rj along 
the plane and d perpendicular to it, where the Josephson current j z is of the 
order of the critical current. In this region, which plays the role of the vortex 
core, the London model fails. Away from the core the streamlines of the 
shielding supercurrents, which also represents contours of constant magnetic 
field, are elliptical except for the zigzags due to the intervening insulating 
layers(Fig.2). It is not crucial what kind of transition is inherent in the mag- 
netic subsystem: spin-flop, metamagnetic or paramagnetic. The main fact is 
the appearance of a magnetic domain with an induced magnetic moment M 
along the vortex. We make the above mentioned conjecture by assuming |J 
that the magnetic moment is constant across the domain 




|M| = 

where p m is the dimensionless radius of the magnetic domain in the coordi- 
nate system of the elliptical cylinder (x = x,y = \ c pcos(p, z = X a bP sirup). 
Then, the set of equations (7) have the following solutions M 



B P = CJo(p) + C 2 K (p) for p<p m , (9) 

BAF = VI T \ K 0(P) fOT P > Pmi 

where B P and Baf denote the magnetic induction inside and outside the 
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domain respectively. The constants C\ and Ci are calculated using (8) and 
the flux quantization condition to obtain 



B T Pmh(pm) ~ D (Pm) I o(Pr. 
1 - h{p m ) 



r , D(pm) K o(Pm) + B T [pmK^pJ - 1] 

° 2 ~~ 1 FT — \ ' y L[J ) 

D (Pm) = o„X X -P0 H T- 



2nX ab X c Ko(p m ) 

where K , K±, Jo, Ji denote the modified Bessel functions whereas Bt = 
M+h Ht- Now we can think of M as of experimental parameter. To calculate 
p m one must first derive the line tension of the vortex and then look for its 
minimum with respect to the domain radius. Standard calculations give 

* . (11) 

It is easy to see that almost all flux quantum is captured in the magnetic 
domain. 



3 ACTIVATION ENERGY 

The first quantitative approach toward intrinsic pinning in layered super- 
conductors was based on the observation that the superconducting order 
parameter should have a periodic spatial variation across the layers fl3| . For 



the present considerations, however, the method of critical nucleus developed 
in 



14] is much more convenient. The activated nucleus consists of a kink- 
antikink excitation, that is, a vortex line segment is thrown to the adjacent 
layer, thereby creating two pancake vortices of opposite "vorticity", as shown 
on Fig. 3. The activation energy can be regarded as the energy barrier for 
intrinsic pinning. Depending on the magnitude of the driving current density 
the process may continue as the single vortex activation or the activation of 
the vortex bundle. First, consider the activation of a segment of a vortex to 
the neighboring interlayer spacing (Fig. 3). The energy associated with this 
process can be written as: 

U = 5E + V k ,-k(R) - U ~ 3o)<PadR. (12) 

The subscript "a" or "b" (of U, 5E and j ) indicates that this quantity is 
calculated for a or b direction in the plane. SE is the amount of condensation 
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and magnetic domain energy that is lost at two points of the layer threaded 
by the kinks separated by a distance R (see Fig. 3). Vk-k(R) is the kink- 
antikink interaction energy. The term proportional to the driving current 
j is due to the Lorentz force. The term proportional to jo is the energy 
associated with the distortion of the line due to the formation of nucleus. 



This term can be estimated from the simple considerations [14 



I f ( du 

joVo d ~ o / dydzC(y, z) 



2 J ~»-~~\»>~>y dz j > 

where the Fourier transform of the compression modulus is given by (15 

B 2 



C{k y , k z 



/i (i + \i b ki + xiki) 



By taking dydz ~ (p /B,u z ~ d, ~ k z ~ k y (X c /X a b) the integral can be 
estimated as follows 

Bd . , 

30a = TT2". (13) 

In the b direction, however, we have an additional contribution from the 
magnetic domain dydz ~ 5ip /8B T according to (11), so we get 

5dB T 

JOb = JOa + 2 • (14) 

iZ ° A ab 

As the current density j drops below jo a single- vortex line can no longer be 
activated due to the confinement energy provided by the vortex lattice. The 
energy 5E is calculated from (7) with the right-hand sides representing the 
vortex cores: 



IN > 77>2/ = 0,z = oj , < ^-,y = Q,z= -d| 
x = ±— ,y = 0,0 < z < —d> . 

The solution is then substituted to the free energy functional (4). Taking the 
limit R — > oo we exclude the energy of the kink-antikink interaction. The 
calculations are involved, so we write down only the results. 

6E a = 2de \n^, (15) 

Safe 
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5E b = de 6 ln-^- , 

Safe 

where e = (^ 2 fi \ 2 ab ), e b = ge In \<p /(irr]B T 



(16) 



The energy of kink-antikink interaction was calculated in |14] 

V K> - K (R) - 



—— -/(— ) 

2A a b A c 



(17) 



where /(t-) 



(\ C /R) - H r j/Lb) for Tj « R « A c 
2( y \ c /Rfexp(-R/X c ) for R»\ c 



We introduce the Ginzburg-Landau critical current density jgl — ^o/i^o^ab^^) 
and the quantity I a;b = 2(j — jo a ,b) / Ugl^V^), then (12) can be rewritten in 
the following way 



Ua 



f, rj T R d „,R S 
2de Q In f- + 1 - f( ) 

I Safe Safe 4A ab A c 



^0,]n-^ + 2 eo / b A_^L /( ^ 

Safe Safe ZA afe A c 



The critical size of the nucleus R c is given as a minimum of (18) with respect 
to R. In the approximation rj«R«\ c corresponding to the current 
regime ( C d/X 2 ab « I a , b « £Jd we get 



n ca,b 



p 2 

S afe 



d 



44,fesV 



2de \ In 



.Cfe, 
r,- x 

, 4^ afe , 



Idla 



d{e b \n[^\-2e 0] j^ 



For the opposite case R>> X c and i c d/\ 2 ab » I a ,t 



(19) 



R 



ca,b 

u c a 



A r In 



'Jtc 
JaK 



2de n ln U 



afe . 



d e 6 ln [±L 

I \Safe, 



Aafe-^a 
Aafe-^fe 



111 



S C c 



hi 



V^aA^, y 
JbKbJ 



(20) 



When the driving current drops below jo the critical nucleus is 3D object 
(Fig.4). In our case it is a parallelepiped of the height R along the bundle 
and of the section S across it. The activation energy is a sum of the volume 
energy due to the Lorentz force and the surface energy. 

U a>h = -jBdRS + 5E a , b (^pj + j 0a , b dR^/BS Vo . (21) 

The second term is the loss of condensation energy (and magnetic domain 
energy in the case of b direction) on both surfaces perpendicular to the 
bundle multiplied by the number of vortices threading these surfaces. The 
third term is the elastic energy released in the surface parallel to the shifted 
vortex joafidR^Q multiplied by the number of shifted vortices \J BS/ip 
(one vortex per plane). The critical nucleus is then S c = ((p / B)(j 0citb / j) 2 , 
R c = 5E a ^ b /(jd(p ), and the activation energy is 

K b = 5E ^{^f) ■ (22) 

4 MOTION OF THE FLUX 

The resistive mechanism in the mixed state is determined by the activation 
process leading to magnetic flux motion (creep). This motion induces electric 
field which can be observed on the current- voltage characteristic. We consider 
the motion of activated kinks along the layers of the length L along the 
magnetic field direction. Assume that each one can reach the boundary of 
the sample before the new one is created. The mean electric field associated 
with this motion is given by 



E = BPLdS c , (23) 

where P is the activation probability per unit volume and unit time. For ther- 
mal activation this probability is given by P = aexp (—Uc/ksT) . There 
is however a crossover temperature T [[16] below which quantum tunneling 
of vortices is dominating. The probability for quantum tunneling is given 
by P = (3 exp (—U c /kBTo) and remains finite even for T = 0. The tem- 
perature To depends on the driving current, critical current, upper critical 
field, and resistivity in the normal state. For our purposes it is worth to 
note that the quantity exp (— [/ c //cbTo) will contribute only for driving cur- 
rents of the order of the critical current fTBfl. The Neel temperature for lay- 



ered antiferromagnetic superconductors varies from hundreds of mK (0.6K 
for ErBa 2 Cu307) to several Kelvin (6.8K for ErNi 2 B 2 C) and therefore both 



S 



mechanisms of activation are present in these compounds. The preexponen- 
tial factors and T cannot be calculated in the framework of thermodynamic 
considerations alone. Fortunately, it was shown in |T7| that the activation 



probability of macroscopic quantum excitations is proportional to j 3 . Thus 
we can assume that P = a j 3 exp (—U c /k B T). Now we can calculate the 
cur rent- volt age characteristics for the current density regimes considered pre- 
viously. For j « jo 



E a = Vo dLa Jo a j exp { -j-^ ( j } , (2 1) 

E b = (p dLa jQ b j exp 

This almost linear dependence of E on j indicates that the resistive mecha- 
nism of bundle activation follows Ohm law. 

For j » Jo and £ c d/\ 2 ab « I a ,b « 




E a = ¥^«oJ expj-— + — W — j, (25) 



P AT 3 J SE b 2de dl h 

2 



For j » j and i c d/X ab » J a6 



= exp |-^ + ^r— ln \w« 1 J ' (26) 

£ fe = (fodLaoj exp <^ -— — + ln 



We can also calculate the rate of flux creep due to the thermal activation 
of vortices. To do this consider hollow cylindrical sample of a radius r and 
the wall thickness I « r placed in the magnetic field B ex > B c \ applied 
parallel to the cylinder axis. The sample has the trapped field E>i n inside the 
hole and trapped flux $ = (B in — B ex )irr 2 . According to the Faraday's law 
electric field due to the change of the trapped flux is equal to (fj, /2)lr(dj/dt). 
Combining this result with (23) we have finally 

BPLdS c + \^lrj t = . (27) 
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This equation can be solved analytically only in the case of the weak 
currents. Consequently for excitations in the form of bundle of vortices (27) 
is written as 



J SEgjb fjo a , b \[ . dj 
^ eXP \-^T [—) ] + dt = °> (28) 

where Q = ^o^oJo / '(Pol) an d 7 = r l/(Ld) is the factor determined by the 
geometry of the sample. The solution of (28) is given in terms of exponential 
integrals and for the case of jo a ,b/j — 1 << 1 it can be approximated as: 

j(0) #(0) k B T (j(0)\ 2 

where 

A(p a 5E ab fj 0a ,b\ J 5E a>b (j 0a ,b\ 2 \ , on , 

= ~^r~ {w)j exv \-^r {w)j I • (30) 

This result is in agreement with the experiments on nonmagnetic HTSC ( see 
||18|| ). For << t << 1/uj the change of trapped flux is linear in time and 
for t » 1/uj logarithmic. In antiferromagnetic superconductors, however, 
we see additional change of characteristic frequency as the magnetic field 
changes its direction in the a-b plane. 

5 CONCLUSIONS 

The current jo plays very important role in all above calculations. For the 
following values d ~ 10~ 9 m, £ c ~ 3 x 10~ 10 m, ^~3x 10~ 9 m, rj ~ 10~ 8 m 
and their typical temperature dependence we can estimate 

^ „ „ IQ-B B I T ° , 

3gl <^ V T c -T 

Since the depairing current is of the order of 10 13 [T c / (T c — T)] 3 ^ 2 then the 
current jo a is of the order of 10 10 i? [T c / (T c — T)] . Although there are no 
precise measurements of spin-flop transition in the antiferromagnetic high 
temperature superconductors, we assume that h q Ht ~ 40mT . The typical 
value of 5.5fi B per rare-earth atom per unit cell gives M ~ 0.37T. It is 
possible now to estimate the change of jo due to the creation of spin-flop 
domain along the vortex : 
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providing that Ht/H c \ — 1 << 1. Similar considerations show that the 
characteristic frequency of flux creep changes of order of magnitude when 
the spin-flop domain is created. 
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FIGURE CAPTIONS 

Fig.l. Schematic drawing of a piece of a layered superconductor. Hatched 
areas (n,n+l) represent superconducting layers with thickness d s . Bold ar- 
rows represent magnetic moments of RE ions lying in the isolating layers of 
thickness dj. The reference frame and the crystallographic axes are shown. 

Fig. 2. Single vortex line lying in the a-b plane and the induced spin-flop 
domain along its core. 

Fig. 3. Single vortex activation. 

Fig. 4. Activated bundle of vortices. 
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